Non-Perturbative Bounds on Hyperfine-Induced Electron Spin Coherence Times 
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We address the decoherence of a localized electron spin in an external magnetic field due to the 
hyperfine interaction with a lattice of nuclear spins. Using a completely non-perturbative method, 
rigorous bounds on the T\ and T2 coherence times for the electron spin are provided. It is shown 
that for magnetic fields B greater than some critical field B c (B c « .001 — 2 Tesla for the systems 
studied here), the z-polarization of the electron spin cannot relax, and hence Ti is infinite. However, 
even at high fields dephasing can still occur. We provide a lower bound on the T2 coherence time 
that explicitly takes into account the effects of a spin echo pulse sequence. 

I. INTRODUCTION 

Isolated electron spin systems are of great interest because of their potential application as coherent quantum 
memory in novel solid state devices[l, 2]. Several quantum computation proposals make use of electron spins, either 
as the qubits themselves or as intrinsic components of one- and two-qubit gates [3-6]. Furthermore, experimental 
efforts are just now succeeding at measuring single electron spin relaxation times [7-9] . Unfortunately, the exponential 
size of the Hilbert space of these spin systems and the variety of channels for decoherence makes an exact analysis of 
spin dynamics virtually impossible (without the use of a quantum computer) [10-13]. As a result, it is often necessary 
to resort to approximations or perturbative treatments to address the timescale of electron spin decoherence, even 
using a simplified Hamiltonian. Because approximations are inevitable, it is important to know which of these can 
be rigorously justified. In particular, it is necessary to know which terms in a very complex system Hamiltonian can 
be ignored because they do not contribute to decoherence on the experimentally relevant timescale. In this paper, 
we consider electron spin decoherence induced by the nuclear spins. In particular, we focus on the relevance of the 
O I S+I- + S-I+ term in the Hamiltonian (see Eq. (1)), which comes from the contact hyperfine coupling of the electron 
and its surrounding nuclei. 

I Previously, the effect of the hyperfine interaction on decoherence has been studied in special cases or using various 
. approximation methods. Recently, Erlingsson et. al used a semi-classical approximation to study electron spin 
correlation functions at low fields (B <C B c , where B c is on the order of the Overhauser field, see below) [14]. They 
found that the correlation functions were characterized by complex dynamics which did not decay, even at long 
times. Work by Khaetskii et al. has characterized the timescale of the logitudinal spin correlation for an electron 
coupled to surrounding nuclei via the hyperfine interaction [10, 15]. In refs. [10, 15], exact solutions for the case 
■ of completely polarized nuclei and for the case of homogeneous coupling, as well as perturbative solutions for the 
general, inhomogeneous situation were presented. From these solutions, the behavior of the general unpolarized, 
inhomogeneous case was inferred. Their results indicate that the logitudinal spin correlation time is on the order 
of HN/A for a single quantum dot and on the order of h^fN j A for an ensemble of dots. Gate imperfections due to 
the hyperfine coupling were also treated using perturbation theory in [16]. However, it was noted in [15] that the 
perturbative solution would diverge to second-order and was therefore potentially unreliable. More recently, Coish 
et al. have treated the behavior of this system using a generalized master equation, again confirming a power-law 
decay for short timescales on the order of HN/A [17]. In general, though, all treatments based on a time-dependent 
perturbation expansion will be valid only at times smaller than some maximum timescale specified by the order of the 
expansion[26]. Although such treatments are useful for examining short-timescale and small-magnitude fluctuations 
in electron polarization, in this paper we are interested in the long-timescale dynamics of the electron spin which leads 
to 1/e decay of the longitudinal or in-plane polarization of the electron. Note that we are not claiming any particular 
functional form for the decay of coherence, exponential or otherwise. Rather, we are evaluating the timescale on 
which the longitudinal and in-plane magnetization decays to a constant fraction (in this case 1/e) of its initial value. 
Moreover, we focus on the high magnetic field regime which is relevant for the recent experimental developments in 
single spin read-out (B = 8 — 14 Tesla in ref. [9]). 

In ref. [11] it was assumed that at sufficiently high magnetic fields the electron-nuclei flip-flop term (S+I- + S-I+ 
in Eq. (1)) which is due to the hyperfine interaction can be neglected as a direct source for electron spin decoherence. 
In this paper, we would like to 1) justify this assumption and 2) determine how strong a magnetic field must be 
applied to safely neglect this term. Rather than resort to approximation to estimate the numerical value of T\ and 
T2, this paper will instead place rigorous, non-perturbative lower bounds on the spin relaxation time of the electron 
in the high-field regime {B > B c ). Furthermore, we will explicitly take into account the application of the spin echo 
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sequence which is often used experimentally to remove inhomogeneous broadening. To determine the importance of 
the electron-nuclear flip-flop term, we use the Hamiltonian presented in [10, 15, 17] which contains only the contact 
hyperfine term. In doing so, we can effectively isolate the decoherence contribution that derives exclusively from 
the hyperfine interaction. When the timescale of this hyperfine-induced decoherence is significantly greater than 
the timescale of decoherence due to other mechanisms (such as spectral diffusion induced by inter-nuclear dipolar 
coupling[irj) then it seems reasonable to ignore the direct effects of hyperfine coupling as a channel for electron spin 
decoherence. 

It is important to define what we mean by "electron spin decoherence". Typically electron spins are assumed 
to undergo simple exponential decay governed by the Bloch equations. Thus the coherence times T\ and Ti can 
be defined as the timescales of this exponential decay. However, when decay is non-exponential, as in this system 
(see [10]), there is some flexibility in the definition of the coherence times. Electron spin dynamics includes both 
short-timcscale fluctuations (which are generally small in magnitude [17]) and long-timcscalc (i.e., 1/e) decay; in this 
paper, we take "electron spin decoherence" to mean the latter. In other words, the coherence time of an electron is 
defined as the time it takes for the longitudinal or transverse polarization of an initially polarized electron spin to 
decay to a specific fraction, e.g. 1/e, of its initial value. 

The decay of the in- plane magnetization is understood to occur in the context of a spin echo experiment. Spin 
echo experiments are designed to remove the effects of inhomogeneous broadening due to variations in the local 
field experienced by an ensemble of electron- nuclei systems (i.e. multiple independent systems). When applied to 
a single electron spin, a spin echo can be understood as removing the inhomogeneous broadening resulting from an 
ensemble of initial states of the single electron-nuclei system. The remaining decay of the in-plane magnetization 
comes from fluctuating Overhauser fields, which cannot in general be removed by the spin echo (but see [18] for 
numerical evidence that a substantial part of hyperfine-induced fluctuations may, in fact, be reversible). Because 
only bounds on relaxation behavior are desired, our analysis of this decay can be completely rigorous, despite the 
complexity of the problem and the size of the Hilbert space. Our results are applicable to localized electron spins in 
semiconductors such as quantum dots and donor impurities. 

This paper is organized as follows: Sec. II outlines the general structure of the problem and presents the results of 
our analysis. The derivation of the theorems used in Sec. II can be found in the Appendix on EPAPS[27]. Applications 
to electron spins in semiconductors of relevance to solid state quantum information processing are made in Sec. III. 
Discussion and conclusions follow in Sec. IV. 



II. RESULTS 



Let S be the electron spin precessing in a magnetic field B applied in the z direction. The electron spin is coupled 
to N nuclear spins Ii, I2, . . . , I at via the hyperfine interaction, characterized by the coupling constants A\, A2, . . . An- 
Due to the hyperfine coupling, the nuclei induce an effective Overhauser field on the electron with magnitude J2j Ajlj. 
Here we consider the Hamiltonian 

H = 1S BS Z + lT B ]T I jz + ]T Aj (s z I jz + \ {S + Ij- + S-I j+ )\ , (1) 

33 

where 7s and 7/ are the gyromagnetic ratios of the electron and nuclei, respectively (7s ~ 10 3 7/). Eq. (1) can be 
written as 

H = H + V, (2) 

Ho = ls BS z + 1I BY,Ijz+Y, A i SJ i^ ( 3 ) 

3 3 

3 

This Hamiltonian conserves the total spin angular momentum J z = S z + 5^ • Ij z ; hence, the total number of "down" 
spins is conserved. For convenience, we will label each of these blocks by the quantum number L, where L is the total 
number of "down" spins (i.e. J z = N — 2L, L = (N — J)/2). Then the Hamiltonian can be block diagonalized based 
on the J z operator [17]. There will be N + 2 blocks, one for each possible value of L, and the L th block will be an 
( N l 1 ) x ( N l 1 ) matrix. For the remainder of the paper, it will be assumed that all analysis applies to some subspace 
labeled by a particular value of L. 

The unperturbed Hamiltonian, H a , commutes with the S z operator. Thus, the eigenstates of H a can be labeled by 
the z-polarization of the electron. On the other hand, the full Hamiltonian H does not commute with S z because V 
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connects electron spin-up states to electron spin-down states, and vice- versa. At high fields greater than some critical 
field B > B c , the Zeeman energy of the electron will dominate the energy of the eigenstates, and the electron spin-up 
and electron spin-down manifolds will be well separated in energy. Thus, if the perturbation is in some sense smaller 
in magnitude than this energy gap, the effects of the perturbation will be small and S z will still be an approximately 
good label for eigenstates of H (see Figure 1). 

Our analysis in Sec. II. A demonstrates that if the magnetic field B is greater than the critical field B c , then the 
electron spin magnetization (S z ) never changes appreciably. Because the eigenstates of H are nearly eigenstates of 
S z at high fields, the z-polarization of the electron does not undergo substantial (i.e. 1/e) decay (this result matches 
the perturbative result found in [15]). Furthermore, based on these arguments, a field-dependent bound is placed 
on the maximum fluctuation of S z from its initial value: the greater the field, the smaller the maximum possible 
fluctuation. Again, it should be emphasized that there may be small magnitude fluctuations in (S z (t)}, which occur 
on a finite timescale. However, because our definition of T\ (the timescale for 1/e decay of (S z )) quantifies the 
large-scale coherence properties of the electron spin, the absence of large-scale decay renders T\ effectively infinite. 
Although Ti is infinite for large magnetic fields, dephasing processes can still occur provided that they conserve the 
z-polarization of the electron spin. Although this statement seems like a contradiction given that the perturbation, 
V, flips the electron spin, virtual transitions involving V 2 can nevertheless flip pairs of nuclear spins without affecting 
the electron spin. For this reason, it is also important to place a lower bound on T 2 . 

Derivation of the lower bound on T 2 in Sec. II. B follows a different methodology than the bound on T\. Because T 2 
is defined as the timescale of the decay of the spin echo envelope [19], a lower bound on the magnitude of (S + ) after 
a pulse echo sequence implies a lower bound on T 2 . Specifically, the standard definition of T 2 is the minimal time, t 
such that, 

|(5 + (2i))| = l (5) 

where a it pulse is applied at time t (we have dropped a factor of h for convenience) . It is assumed that for an infinite 
system, there will be no large-scale recurrence in S+(2t), but only dissipative behavior. Small magnitude recurrence, 
due for instance to echo modulation [20], may still exist, but large-scale recurrence on the order of 1/e is assumed to 
be absent. Even if this large-scale recurrence does exist, the lower bound on T 2 is still valid. In particular, if it can 
be shown that for some value of t', \(S+(2t'))\ > l/2e, it follows that T 2 > t' . 

This bound on \(S + (2t')}\ can be proven using a variational argument. Every trial state \(j)) induces a probability 
distribution over the eigenenergies of H. The mean of this energy distribution is given by the expectation value 

(E) = (<p\H\<f>). (6) 

Likewise, the variance of this energy distribution is 

((AEf) = (4>\H 2 \4>)-^\H\ ( t > ) 2 . (7) 

The variance of the energy will be equal to zero if and only if \(p} is an eigenstate of H (or a superposition of degenerate 
eigenstates). For all other states which are superpositions of non-degenerate eigenstates, the variance will be greater 
than zero. The essential feature of the theorems given in Sec. II. B is that the dynamics undergone by any state \(j>) 
are limited in timescale by the variance in energy of that state. Using this idea, it can be shown that the dephasing 
processes of the system have a timescale governed by the energy 

Conceptually, our bound on T 2 makes mathematically rigorous the idea that the coherence time of a state is inversely 
proportional to its linewidth. 

Before proceeding to the proof of these results, two important points must be made. First, this method constructs 
only a rigorous lower bound on T 2 ; it does not predict the actual value of T 2 nor does it claim to be a tight lower 
bound on T 2 . In fact, numerical work suggests that the spin echo pulse may actually remove nearly all hypcrfinc- 
induced decoherence, leading to only small magnitude decay of S+(t) under the spin echo sequence for all times, t [18]. 
However, this evidence is numerical and has been confirmed only for small numbers of spins (N ~ 10). In this paper 
we are interested in what rigorous statements can be made regarding the T 2 time for a large system (N ~ 10 4 — 10 6 ). 

Second, it is the special structure of this system at high B fields that allows a useful bound to be constructed. 
Though this method is general and can be used to construct valid bounds on T 2 for other systems, these bounds 
may prove to be useless in practice (for instance, the rigorous bound that T 2 > 10~ 30 s is valid, but is a practically 
useless bound). Fortunately, we find that our method provides practically useful lower bound values for real systems 
of relevance to quantum information processing (see Sec. III). 
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A. Bound on T\ 



The proof of the bound on 7\ will be given first and will be organized as follows. First, we will show that the 
eigenstates of H can be divided into two subspaces, the + subspace and the — subspace (Theorem 1). Eigenstates 
contained in the + subspace contain large contributions from electron spin-up states, and eigenstates in the — subspace 
contain large contributions from electron spin-down states. In other words, S z is an approximate quantum number 
for the true eigenstates of H . Second, we will show that this fact places bounds on the eigenvalues of H (Theorem 
2). Next, we will show that every eigenstate of Hq is almost completely contained within cither the + or — subspace 
of H (Theorem 3). To complete the bound on Ti, we will show that these previous statements imply a lower bound 
on the decay of S z from its initial state (Theorem 4). 

For future convenience, we define 



A[k] 


= E ^ 


(9) 




j£k largest 




A 2 [k] 


= E ^ 


(10) 




j£k largest 




A 


= E^ 

3 


(11) 


A2 J max 


= E^ 2 - 


(12) 



In other words, A[k] is the sum of the k largest hyperfine coupling constants. A max is the sum of all N hypcrfinc 
coupling constants and therefore represents the maximum magnitude of the Overhauser field. The matrix Hq com- 
mutes with the S z and Ij Z operators; hence eigenstates of H can be labeled by the z-polarization of the electron 
spin and the z-polarization of the nuclei. Let Iff) and |JJ.) represent the +z and — z polarized states of the electron 
spin, respectively, and let z be an A^-bit string of +ls and —Is representing the z-polarization of the N nuclei. The 
eigenstates of H are of the form Iff, z) or |JJ-,z') (it should be remembered that z will contain L down spins and z' 
will contain L — 1 down spins, since the total number of down spins, including the electron spin, is conserved). The 
eigenvalue Eq z corresponding to the eigenstate |ff, z) is 

E L = +\^lsB + \(N- 2L) h 7l B + \h 2 ]T Aj Zj . (13) 

3 

Similarly, the eigenvalue Eq , corresponding to the eigenstate |JJ-, z') is, 

<z< - -\nisB + \{N-2L + 2) h^B \h 2 ]T Ajz'j. (14) 

3 

From Eqs. (13) and (14) it is clear that the eigenvalues of H fall within certain ranges. In particular, we can write 

E^ < E^ < E^ 



EX ■ < Ex , < E n 

0,min — 0,z' — \),max 



where we have defined 



0,7m ro 


= +±h ls B + ^(N- 


2L) h^B - X -fi 2 (2A[L] - A max ) 


(16) 


O.max 


= +±h ls B + ^(N- 


2L) ft 7/ 5 + ^H 2 (2A[N - L] 




(17) 


F 4 

O.mm 


= - l -h ls B + l -(N - 


2L + 2)h 1I B-^Ti 2 (2A{N- 


-L + l}- A rnax ) 


(18) 


O.max 


= -±h ls B + \{N- 


2L + 2)Ti- fI B+hl 2 (2A[L- 


l] -^-max) • 


(19) 



It is critical to all proofs that follow that there is some energy gap AE between the electron spin-up and spin-down 
subspaces (see Figure 1). In other words, the lowest energy electron spin-up state must be greater in energy than the 
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highest energy electron-spin down state. Using Eqs. (16) and (19), it can be demonstrated that an energy gap exists 
if 



His - li)B > h 2 (A[N - L] + A[N - L + 1] - A max ) . 



(20) 



An additional condition that emerges from Theorem 2 sets the critical field slightly higher than Eq. (20) would imply. 
Furthermore, we would like the definition of B c to be independent of L. Therefore, for the remainder of the paper, 
we will assume that we are in the high-field regime defined by 



B > B C1 



where 



B r 



hA r , 



Is - li 



(21) 



(22) 



Typical values of this critical field for various systems are given in Table 1. 

Using these facts, we now show that the eigenstates of H are near-eigenstates of S z . We define the projection 
operators and such that 



where Xj is the identity operator on the nuclear spins. Then the following theorem can be proven: 

Theorem 1. The eigenstates of H can be divided into two subspaces, which we will label + and 
in the + subspace obeys the conditions 

(V + |n t |^+) > Ci 

and an eigenstate \tl>~) in the — subspace obeys the conditions 



(r 



IT 



where 



U,2 



1 1 

2*2, 



4U 2 



Jj. ^ 

O.max 



(23) 
(24) 

An eigenstate 

(25) 
(26) 

(27) 
(28) 

(29) 



such that 1 — > +, 2 



It is useful to note that 



l -h 2 ^A[L]A[N- 



l-L]. 



(30) 



(31) 



for all L. Above the critical field given in Eq. (22), we find that £i > 1/2 and £ 2 < 1/2- Hence, we can write the 
eigenstates of H as either \4>t) or where the + eigenstates are those which contain contributions primarily 

from electron spin- up states (Iff)), and the — eigenstates contain contributions primarily from the electron spin down 
states (IU-)), 



<^ + |n t K + ) > 1/2 
(^"|n^|V fe -> > 1/2. 



(32) 
(33) 



To summarize, we have shown that at high fields, states can still be approximately labeled according to their electron 
spin polarization. 

We now proceed with the second step of the proof. 
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Theorem 2. The eigenvalue corresponding to eigenstate an< ^ ^ e eigenvalue E k corresponding to eigenstate 
| V'fe ) obey the following conditions 

E+ nin < E+ < E+ ax (34) 

and 

E mm < E- < E max (35) 

where 

E>min = -^0,min — V ma x (36) 

Emax = ^max + Vmax (37) 

E m in = ^0,min — Kraax (38) 

E m ax = -^0,max + Vmax- (39) 

Theorem 2 demonstrates that the eigenvalues of the + and — eigenstates of H fall in approximately the same ranges 
as the eigenvalues of the ft and -I)- eigenstates of Ho (c.f. Eqs. (16-19)). 

The third step in the proof is to show that any eigenstate of S z , in other words, any state with a well-defined 
electron polarization, is contained almost completely in either the + or — subspace (note that this statement does 
not immediately follow from Theorem 1, in which we proved the converse; namely that + or — eigenstates of H are 
almost completely contained within the ft or JJ. subspace.) 

Theorem 3. Let |ft,^>) be an arbitrary state in which the electron is polarized in the +z direction, i.e. \ip) is some 
arbitrary state of the nuclear spins. Then, 

(ft,V|n_|ft,v) <4, (40) 

where II is the projector onto the — eigenstates of H and 



e. 



-Bjf +v 



2 x 0,max 0,min J * max 

J^max a Q,min I 



(41) 



Similarly, 



(^|n+|^><4, (42) 

where II_|_ is the projector onto the + eigenstates of H and 

( e4 _ c4 A 2 I T/2 

„ \ O.mai O.min ' y max 

4 = — 2 — 43 

1 min O.max J 

To simplify future analysis, we will use 

e = max{e t ,e 4 } (44) 

in place of and ejy. Then Theorem 3 demonstrates that an arbitrary state with a ±z polarized electron has at most 
an amplitude of e on the =F subspace. 

Finally, let us consider the time evolution of the state |ft, tp). We will show that at high fields, the z-polarization 
of a state docs not vary substantially from its initial polarization. In other words, at fields B > B c , T\ is infinite 
because the ^-polarization of the electron does not undergo substantial decay. 

Theorem 4. Let U be the time evolution operator given by, 



U = exp(iffi). 



(45) 
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Given an arbitrary initial state with the electron polarized in the +z direction, the following inequality holds for 

all times t, 

{it, 4>\ U\t)S z U{t) % V) > \ (1 - 8e 2 ) . (46) 
Similarly, given an arbitrary initial state |JJ., ip) with the electron polarized in the —z direction, 

% i,\ UHt)S z U(t) |4, V) < ~\ (1 - 8e 2 ) . (47) 

Given that the expressions for e 2 and £2 are complicated, it is useful to observe their general form. If we assume 
that h~j s B > h A 

max, we can expand e and £2 to lowest order in 1/B. To the lowest order in 1/B we find that 

e 2 w 5C 2 /4 and 



& » 777^ ^-2 (48) 



x -h A A 2 

e 2 and £ 2 are essentially the squares of the ratio of the Overhauser energy h~ 2 A max to the electron Zeeman energy 
hjs'B. If this ratio is small, then the z spin of the electron undergoes only small fluctuations on the order of e 2 . Thus, 
for fields greater than B c , decay is suppressed and 7\ is consequently infinite. 



B. Bound on T 2 



Let us now turn our attention to the second relaxation time T 2 . T 2 describes the decay of the in-plane magnetization 
under a spin echo experiment, i.e. the decay of the expectation value \(S + (2t))\, where a single spin echo 7r-pulsc is 
applied at time t. The bound on T2 is based on the variational principle mentioned above. It will be shown that 
the + and — cigcnstates obey very similar evolution equations. Thus, the dcphasing which occurs between + and — 
eigenstates is limited by the "magnitude" of this difference in dynamics. 

We first begin with some general observations. Let \<p) be a wavepacket formed from a linear combination of 
eigenstates of the Hamiltonian H . Then we can write, 



(E) 


= (4>\h\4>) 


(49) 




k 


(50) 




= ^2p k E k 

k 


(51) 


(E 2 ) 


= (4>\H 2 m 


(52) 




k 


(53) 




= Y,p* E l 


(54) 



k 



where \^k) is an eigcnkct of H with eigenvalue E k and where we have defined p k = |(i/'A;|0)| 2 . We can view p k as 
a probability distribution over the eigenstates of H induced by the state \<j>). Then we can describe the root mean 
squared width AE of the distribution of energies as, 

AE 2 = (E 2 )-(E) 2 (55) 

= J2Pk( E k-(E)f- (56) 

k 

Theorems 5 and 6 will apply this idea to bounding the dynamics of \<j>). 

Theorem 5. Let \<f>) be a trial state with energy width AE, let AE max — aAE be some energy width cutoff with 
a > 1, and let M be the set of eigenstates of H states such that \E k — (E)\ > AE max . Then, 



(57) 
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Theorem 6. Given a state <fi with energy width AE, then for times t which satisfy, 



Kir 

1 < 2Z*E> ™ 



with a > 1, the free evolution of \(f>) is given by, 



U(t) |0) = exp (iO)y/l - A 2 10) + A |r) , (59) 

where 9 is some arbitrary phase, \r) is a residual vector orthogonal to \<j>) and the parameter A satisfies the inequality, 

A 2 <c(AE,t) 2 . (60) 

The parameter c(AE, t) is a function of the root mean squared energy width AE and the time t such that, 

, A „ / 1 \ { aAEt\ 1 , . 

c(Ai^ = l-(l-^c s(— J + ^, (61) 

and a is a parameter that can be chosen to optimize the resulting bound on T 2 ■ 

Using Theorem 6, it is possible to bound the dynamics of any initial state \<j>) by evaluating its energy width, AE. 
We will come back to the choice of a later when we evaluate the T 2 bound for specific systems in Sec. III. 

We now turn to the problem of using the results of Theorems 5 and 6 to place bounds on T 2 for our system of an 
electron spin coupled to a lattice of nuclear spins. We assume that we start with a thermally equilibrated ensemble 
of nuclear spins. A spin echo experiment is begins with an initial it/ 2 pulse denoted by the operator, 

fl*/2 = (KX1H + H=><W)®z} (62) 

where 

k> = ^m-m (63) 

m = ^m + m). (64) 

The system then undergoes free evolution for a time t. This free evolution is followed by a tt pulse, denoted by the 
operator 

R« = m(n\ + \n)m®ii, (65) 

which flips the electron spin. Finally, there is a second period of free evolution for a time t. This sequence is described 
by the unitary spin echo evolution operator, U, 

U = U(t)R n U(t). (66) 

The decay function v(2t) is given by the magnitude of the in-plane magnetization after this sequence. 

At high fields, it is possible to cool the system well below the Zeeman energy of the electron, thereby approximately 
polarizing the electron in the — z direction. However, because the Zeeman energy of a nucleus is 10 3 smaller than 
that of the electron, the nuclei will, in general, be highly thermalized [11]. For the purposes of our analysis, we can 
then assume that in the initial state of the system, the electron is polarized in the — z direction (i.e. aligned with the 
field) and is uncoupled from the nuclei, which arc in a completely mixed state. The density matrix for such a system 
is given by 

Po = fc|^)(^|®2:/ (67) 
where k is some normalization constant which we will drop for convenience. Thus, we need to calculate the quantity 

v(2t)= Tr\Rl /2 U^S + UR w/2PQ \ . (68) 
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Using the cyclic property of the trace and Eq. (67) , we can simplify this expression to 



v(2t) = 



Tr[mS + UR w/2 p Rl /2 
Tr\iflS + U(Is + R,)®Ii 



Tr 



Tr 



Ws+u 



+ 2 Tr 



(69) 
(70) 
(71) 
(72) 



where Is is the identity operator on the electron spin. 

If we evaluate the trace in Eq. (72) using the eigenbasis of H, we can separate the contributions of the + and 
eigenstates, 



w(2t) = ^Tr^W S+URk 



> 



J2(^\&S+UR V \^) 



Y,{^t\U ] S + UR^\^i) 



(73) 
(74) 

(75) 



Using Theorem 1 and Theorem 4, we find, after some tedious algebra, that the contribution of the + eigenstates is 
negligible, 



^(il>+\WS + UR«\4>f) 



< 2eV6+4e 2 +6- 



(76) 



If we let | tp i ) 
— states, 



fr.i 



+ 



JJ-, ipjy i ) be an eigenstate in the — subspace, we can also simplify the contribution of the 



> 



E(ft>V^ &S+UR* 



2eV6-4e^2- 



Combining Eqs. (75-76) we obtain 

1 



v(2t) > 



-2eV6-2e 2 -i6-2e 2 6 



The problem of evaluating v(2t) has been reduced to determining the time evolution of the states 



(77) 



(78) 



and 



"fr' V'jj. ij- The physical interpretation of the derivation of Eq. (78) is that we have separated out terms which primarily 
contribute to dcphasing processes (i.e., the terms in the sum over i) from terms due to the longitudinal relaxation 
of the electron (i.e., the remaining terms). These longitudinal terms are on the order of £2 or e 2 since this is the 
magnitude of longitudinal relaxation demonstrated in Theorem 4. 

We must next evaluate the sum in Eq. (78). It is clear that under free evolution, 



JJ-, f/'jj. ij simply acquires an overall 
phase (minus some small residual), since it is the projection of an eigenstate of H onto the JJ. subspace. This statement 



is formalized in Theorem 7. A similar idea can be used to characterize the evolution of 



. Because 



is 



not necessarily close to an eigenstate of H, its evolution can be complicated and may not in general be characterized 



simply by an overall phase. If we could calculate the energy width of 
its dynamics. However, because the explicit form of the nuclear state 



ft, ip„ i ), then we could use Theorem 6 to bound 



1>. 



Jj-. t 



is not known, we must use a more subtle 



method of obtaining the energy width. Theorem 8 shows that a trial state \<j>) can be constructed with energy width 
given by a quantity ||TU|| 2 that is determined by ^2, mas, E^ ax , and E { 



4 

Q,max ' 



It is also shown that 



projection of \<j>) onto the ft subspace. Using this fact, Theorem 9 then bounds the evolution of 
present these three results. 



1t,ip u ) is the 
ft>^,i/- Wc now 
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Theorem 7. Let ) = ft,^ j) + ■0-)^ 4 ) ^ e ari eigenstate of H in the — subspace. Then 



U 



r i\.,i 



where 



Again, because 



ft,^. 



1t, i 



(79) 



(80) 



J|, ipy i ) is close to an eigenstate of H, its evolution will essentially be characterized by an overall 



phase plus a small residual. Because 



ft, / is not necessarily close to an eigenstate of 77, we cannot use the same 



reasoning to characterize its evolution. Instead, we show that 
we can characterize. 



ft, V'ju) is close to a trial state \<j>i) whose evolution 



Theorem 8. There exists a trial state \<pi) such that 

n t l^i) = |ft,^) 

and 



(^|iJ 2 |^)-(^|7f|^) 2 < \\W\\t 



where 



\W\\ 2 < 



h 4 



4 (^E m ax + -2'0,max) 



.max (l + e) 



(81) 



(82) 



(83) 



According to Theorem 8, the trial state \<pi) has an energy width of ||IU|| 2 . Using Theorem 6, we know that for 
t < hTr/2\\W\\ 2 the evolution of \<pi) will then be well-described by phase evolution with some arbitrary phase denoted 

by (corresponding to 9 in Theorem 4) plus some small residual. Because \<j>i) is close to the state ft, V^,^) 



Theorem 9 now uses the same reasoning as in Theorem 7 to characterize the evolution of 



Theorem 9. Lei \tp i ) 



ft,^. 



JJ-,'011. A be an eigenstate of H in the — subspace. Then for t < K / k/2\\W\\ 2 , 



U 



w/iere £/ie following inequalities are satisfied, 



\ 2 t <c(\\W\\ 2 ,t)\ 



(84) 



(85) 



(ii-,h,i\^h,i) <£ 2 



< e. 



(86) 



(87) 



Using Theorems 7 and 9, we obtain the following results for the evolution of 
echo experiment: 



l$-,if> hi ) and 



U 



under the spin 



(88) 
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u 



= e i K, I +"t, l )*y / l - A? 
+ (e^SUR^ - UR 



One effect of the spin echo evolution is especially important to note. If we compare the free evolution of 



(89) 



and 



1h V^.i/i we nn< i that they acquire phases of cJ^^t and w^i, respectively. Thus, there is a relative phase difference 
between the two evolutions, which would lead to a fast dephasing time T 2 when we sum over i in Eq. (78) to account 
for our initial, incoherent thermal distribution. However, when the spin echo experiment is performed, Eqs. (88) and 
(89) show that both states evolve with the phase factor of (w^ i + w^Ji, leading to no relative phase difference. As 
expected, the inhomogeneous dephasing due to an incoherent distribution over initial states is removed by the spin 
echo experiment. 

Finally, we can use Eqs. (88) and (89) to obtain, 



ws+u 



ti^li) >l-h-ec(\\W\\ 2 ,t)-2c(\\W\\ 2 ,ty. 
Combining Eq. (78) and Eq. (90), we obtain as a final result 

v(2t) > I (l - 26 - 4e 2 - 4eV6 - - e c(\\W\\ 2 ,t) - 2c(\\W\\ 2 , tf) 



(90) 



(91) 



for t < hir/2a\\W\\ 2 . 

Though this bound is rigorous, it is somewhat complicated. In practice, numerical evaluation of the expressions for 
£,£2, and ||W|| 2 can be greatly simplified if we assume that B ^> B c , as in Eq. (48). In this case, we can also use the 
approximation that 



\\W\\ 2 *W„ 



where 



W„ 



h 4 A 



2,max 



A{h ls B) ' 



(92) 



(93) 



It should be emphasized that this approximation (B 3> B c ) is not intrinsic to our proof and is used only to make our 
results more compact and easier to present. If a strict bound were desired, we could drop this approximation and 
keep only the weaker condition that B > B c . 
To obtain a bound on T 2 , we set v(2t) = l/2e, 



1 



1 



c(||W|| 2) tf + ec(\\W\\ 2 ,t) + 6 + 26 2 + 2 ev /£ 2 + 2e% + - - - = 



2e 2 

Let c max be the positive root of this equation. Then we solve for % in terms of a such that c(||W|| 2 ,i6) = c„ 



U = max 

a 

This value of tb will be a lower bound on T 2 . 



AEa 



cos 



1 



(94) 



(95) 



III. APPLICATION TO PHYSICAL SYSTEMS 



To demonstrate how our bounding method can be used in practice, we apply it here to three systems of relevance 
to quantum information processing[3, 21]: (1) a donor impurity in Si, (2) a donor impurity in GaAs, and (3) a GaAs 
quantum dot. Table 1 shows the important physical parameters for these systems. To obtain a lower bound on T 2 
for each system, we must evaluate Eq. (95) for the given system's parameters. Because we are free to choose the 
parameter a in Theorem 6 with only the restriction that a > 1, we can select a to maximize tb- In each of the 
following examples, this optimization procedure was carried out to obtain the largest possible tb- 
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Consider first the case of a phosphorous impurity atom in a lattice of natural silicon (which contains 4.67% of the 
spin-1/2 29 Si isotope). Spin coherence is affected by hypcrfinc coupling between the electron and nuclei and dipolar 
coupling between nuclei. The hyperfine coupling constants are given by the equation [11, 22], 



where Xj is the position of nuclear spin j and W is a hydrogenic orbital for the impurity electron. Using our bounding 
technique, we can determine when the S+i_ + S-I+ terms in the Hamiltonian can be safely neglected as the primary 
source for electron spin decoherence. The coherence time due to dipole-dipole coupling of nuclei for a sample of natural 
Si in the presence of an external field aligned with the [111] direction has been estimated to be T 2 ~ 0.65 msec [11]. 
This value is essentially independent of the magnetic field strength. In contrast, Figure 2 shows that our lower bound 
is directly proportional to B. Using the data from Figure 2, we find that at approximately B = 3.0 Tesla, decoherence 
due to the dipolar nuclear coupling dominates that due to the hyperfine interaction. Thus, for this system, at fields 
above 3.0 Tesla it is safe to neglect the S+I- + S-I+ terms in the Hamiltonian as sources for decoherence. 

Next, we investigate the case of a donor impurity in GaAs and a GaAs quantum dot (20 nm radius, 10 nm 
thickness). Ga and As nuclei are both spin 3/2. We can easily extend our spin-1/2 formalism to include higher 
spin nuclei if we note that a spin-3/2 particle can be treated as a composite system of three spin-1/2 particles. 
However, we also note that the Ga and As nuclei will have different gyromagnctic ratios 771-Ga = 8.16 x 10 7 (s Tesla), 
= 6.42 x 10 7 (s Tesla), and l 75 As — 4.58 x 10 7 (s Tesla). As a result the Zeeman energy of the electron can be 
exchanged with the nuclear Zeeman energy, providing a new channel for T\ relaxation. In order to apply our bound to 
the GaAs systems, we therefore must make the assumption that at high fields nuclear polarization is not transferred 
between nuclei of different types. We then let 77 = y^Gai and apply our proof as before. Again, the important 
parameters are shown in Table 1. Figure 3 shows the lower bound on T 2 as a function of B for the donor impurity in a 
GaAs lattice (solid line) and a GaAs quantum dot (dot-dashed line). A lower bound for T 2 is provided for fields greater 
than ~ 7.5 Tesla. For both systems, our method also shows that T\ is infinite above a critical field of ~ 2.25 Tesla. 
Previous studies have shown that nuclear spectral diffusion causes decoherence on a timescale of ~ 10 /xsec for both 
these GaAs systems [11, 23]. Hence, Figure 3 shows that for magnetic fields greater than approximately 12.5 Tesla, 
nuclear spectral diffusion will be the dominant decoherence mechanism and the S+I- +S-I+ terms in the Hamiltonian 
can be neglected. 



It is useful to assign a physical significance to the quantities that appear in our bounds. For instance, the relevant 
quantity in the T\ bound is given alternatively by e 2 or £ 2 , both of which quantify the mixing of the spin- up and 
spin-down eigenstates of Ho. These quantities scale as the square of the ratio of the Overhauser energy H A max and 
the electron Zeeman splitting, hjsB. In the language of perturbation theory the large energy denominator suppresses 
mixing of spin-up and spin-down states via the perturbation V. 

On the other hand, the operator which sets the scale for T 2 is W, which can be written as, 



This expression motivates us to identify the operator W with difference in effective Hamiltonians felt by electron 
spin-up and electron spin-down states. For instance, the second term in this equation describes the Hamiltonian 
experienced by an electron spin-up state scattered via two applications of V into a second electron spin-up state. On 
the other hand, the first term applies the same scattering process to an electron spin-down state. Taking the difference 
of these two terms allows cancellation to occur, reducing the magnitude of W, and corresponding to the refocusing 
of spin states accomplished by the spin echo experiment. If no spin echo were applied, then the first of these terms 
would be absent, and the magnitude of W would be substantially larger, leading to a faster decay of coherence. 

In summary we have used a completely non-perturbativc approach to place a rigorous lower bounds on the T\ 
and T 2 coherence times of an electron spin coupled via the contact hyperfine interaction to a lattice of nuclear spins. 
For 7\, we have shown that for B > B c , the z-polarization of the electron is nearly conserved by time evolution; 
hence T\ is infinite. For T 2 , we have obtained an analytic expression for a lower bound tb, having explicitly taken 
into account the effects of a spin echo experiment, which would be required to remove any additional inhomogencous 
broadening. This analysis holds above some critical field satisfying the high field condition, B > B c . We have also 
assumed that at the start of an experiment, the electron is polarized and uncoupled from the completely mixed 
nuclear initial state; similar results can be derived for any highly mixed nuclear initial state. In general, this method 



Aj = — ls nm(^)\\ 



(96) 



IV. DISCUSSION AND CONCLUSIONS 




(97) 
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System 


B c (Tesla) 


ft/Wma^s/Tesla) 


Si:donor 


9.2 x 10~ 4 


1.31 x 1(T 3 


GaAs:donor 


2.24 


1.45 x 10~ 5 


GaAs dot 


2.26 


9.67 x l(T fi 



TABLE I: This table lists the relevant parameters for our calculations of the lower bounds on 7~i and T 2 for the electron spin. 
B c is given by Eq. (22) and W ma x by Eq. (93). The parameters for the Si:donor and GaAs:donor systems were calculated using 
simple hydrogenic orbitals as described in [11, 22]. The GaAs quantum dot considered had a radius of 20 nm and a thickness 
of 10 nm. The parameters for the quantum dot were calculated as described in [11]. Here TiB /W max provides an approximate 
bound for T2 for any field strength B. 

can be used to determine when the contact hyperfine interaction can be safely neglected as a direct mechanism for 
electron spin decoherence. We have demonstrated the utility of the bounds for a donor impurity in Si, and -within 
the approximation that spontaneous nuclear polarization transfer does not occur- for a donor impurity in GaAs, and 
a GaAs quantum dot. 

Future directions for this investigation are two-fold. First, numerical simulations indicate that the coherence time 
T 2 due to hyperfine coupling only might, in fact, be substantially higher than the lower bound given here, possibly even 
infinite [18]. The primary contribution to T 2 comes from the diagonal matrix elements of W, which corresponds to the 
self-energy of the perturbation. Removing this term would greatly improve the lower bound on T 2 ; unfortunately, it 
is unclear how this term can be suppressed. Also, we have not taken advantage of the substantial symmetries present 
in the system [18, 24]. Taking advantage of these symmetries in our analysis could dramatically increase our lower 
bound on T 2 . Second, this method might be extended to yield bounds on the fidelity of single qubit operations on 
the electron spin, which would be important for any solid state quantum computing applications. Finally, it is likely 
that our method could bc extended to any quasi-two level system with a significant energy gap that is coupled to 
a bath via an off-diagonal perturbation. It would be useful to generalize the method so that it could be applied to 
other systems of interest to quantum computation. 
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FIG. 1: A schematic representation of the eigenspectra of Ho and H . The eigenstates of the unperturbed Hamiltonian can be 
grouped into two subspaces, ff and JJ. based on the z-polarization of the electron. In Theorems 1 and 2, we demonstrate that the 
eigenstates of the full Hamiltonian H can be grouped into two subspaces, + and — , based on the approximate z-polarization 
of the electron. Theorem 3 further states that the + and — eigenstates fall within the same range as the -ft and JJ. eigenstates, 
respectively. 
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V. APPENDIX 

Proof of Theorem 1. Let \ip) be an eigenstate of H such that, 

H\4,)=E\1>). (98) 

Without loss of generality, we can write \tp) as, 

W = |1h^> + I^>- (99) 

where 

= n t |^), (ioo) 
|^,^) = n^), (101) 

Recall that 11^,11^ are the electron spin projection operators defined in Eqs. (23) and (24). Multiplying Eq. (98) on 
the left by either (ff, V'tl or (JJ., we obtain the two equations, 

(1\,^\H \H,^) + (1\,^\V\^^) = E(it,^\it,^) (102) 
(1$,, ^| H % ^> + (4, ^| V ^) = E (4, V^l 4, H) ■ (103) 

It is straightforward to set a bound on the quantities (-ft - , V'-frl ID", V'-ft-) an( i (-U-i V'JjJ -^o 14, V'Jj-) since the eigenspectrum 
of H is known (see Eqs. (13) and (14)). Specifically, it can be shown that, 

(ti> t \Ho\t^t) >Kmin4 (104) 

and 

{^\H %^)<El max c\, (105) 

where 



Ctf = y/(t,iPit\ It^lt) (106) 

^ = Ji^Mi^h)- (107) 

Here, we pause to introduce several relations involving matrix norms which we will make use of in what follows. Let 
X be an arbitrary Hermitian matrix. Then following the standard definitions, we let H-X"^ represent the maximum 
absolute column sum norm 

II^Hj = maxV \Xij\. (108) 

The maximum absolute row sum norm is given by 

||X|| oo =max^|X ij |. (109) 

3 

Finally, the spectral norm, can be defined as 

|LY|| 2 = max 1(^1 X\cj>)\. (110) 

\<t>) 

When applied to a vector HIV")!^ is simply the L 2 norm of We will also use two useful relations regarding these 
norms [25] 

|LY|V>|| 2 < ||X|| 2 |||V>I| 2 (HI) 

ll*ll 2 < ll^lllll^lloo- (H2) 

We can use these relations applied to the matrix V to bound the remaining quantities in Eqs. (106) and (107). 
Using Eq. (Ill) we obtain 

imk^j v \u.$>n)\ < c. A c„iin A yn„iL. ni.n 
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Because the absolute column sum and absolute row sum norms of V are easy to calculate, we can use them to bound 
the spectral norm of V within a given L subspace using Eq. (112) 

||n^Vn 4 || 2 < h 2 ^A[L]A[N - L+l] (114) 
= V max (115) 

Similarly, 

lin^VTLfrHa <V max . (116) 

Then we obtain the two bounding inequalities, 

E tmin4 - C t C ^ V ma X < Ee\ (117) 
EtmaA + C t C^V max > E C \ ( 118 ) 

(119) 



Combining Eqs. (117) and (118), we obtain the inequality 

c t + H + c t c ^ v 1 >0- (120) 

* max 



Since the normalization of \ip) implies that + c| = 1, some brief algebra then yields 

V O.min O.max 



v 2 

4 2 , v max \ n /101 \ 

c t - c t + T^t — — « ^ °- ( 121 ) 



Solving the equality for c|, we obtain the two solutions = £1, £2 where £1 and £2 are given in Eq. (29). Thus, using 
Eqs. (100) and (106), we have shown that every eigenstate \ipf-) of H satisfies one of the two inequalities, 

(^|n t |Vfc> (122) 

or 

(^|n t |v> fe > <6 (123) 

QED □ 

Proof 0/ Theorem 2. We will show only one of these bounds as the rest can be demonstrated by identical procedures. 
For a + eigenstate, we know from the proof of Theorem 1 that 



H < V6 (124) 
H > V6- (125) 

Then Eq. (117) gives us, 

(126) 





> 


O.min 


\ v max 

V ^ 




> 


E 11 

O.min 


— Vmax 






E+. . 

min 





(127) 
(128) 

Where we have recalled that our high field condition ensures that £1 > £2 (see Eq. (29)). The other bounds are derived 
in exactly the same way. QED □ 

Proof of Theorem 3. We will demonstrate this theorem for the spin-up case; the spin-down case follows the same 
argument. Let IfT,?/') be an arbitrary state with the electron polarized in the +z direction. Since Iff,?/') nas n0 spin- 
down component, its expectation value over V is 0. Given the eigenspectrum of H , it is straightforward to show 
that, 



> (-tt.MHlit.ib) > Et 



(129) 
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Given the matrix norm of the perturbation V in Eq. (116), it is equally straightforward to show that, 

(it, M (H (H)f \t, i>)<\ (El max El mm ) 2 + Vl ax . (130) 
However, we can also write the left hand side of Eq. (130) as, 

(tiP\(H-(H}) 2 \tiP) = Y,{\(^^)\ 2 ( E k - W) 2 + \(tM^)\^E,: - (H)) 2 ) (131) 

fe 

= E (Pt ( E i Wf+Pk (E- k (H)) 2 ), (132) 
fe 

where 

pt = \(itMk)\ 2 - ( 133 ) 

Then using Eqs. (34) and (35), we can bound this quantity as, 

Y / pi(E+-(H)) 2 +p-(E--(H)) 2 > (0) 2 ^p++(E max -El mm ) 2 J2Pk (134) 

k k k 

= (s max -<„ m ) 2 (t^|n_|tv)- (135) 

Combining Eqs. (130) and (135), we obtain, 

(1hV|n_hh^><4 (136) 



where 



I ( _ \ 2 , T/2 

4 I 0,max 0,min ) ~ "m< 



(137) 



(E- -Eft V 

( Cmra ^O.min J 

QED □ 

Proof of Theorem 4- We will demonstrate the inequality for the initial state |fT, ^i>); the spin-down inequality can be 
proved in a similar manner. Given the initial state \if, ip), the state at time t will be given by, 

m))=U%1>). (138) 
The amplitude in the electron spin-down subspace after time t is given by, 

||IWft,V>|| 2 = max|<^|£/|ft,V>| (139) 
m 

= max|(^|(n + + n_)[/(n + + n_)|tV'>l (140) 
\<p) 

< max <j>\ n+C/n+ \if,^)\ + max |<^0|n_tfII_ (141) 

\<t>) 10) 

< max||n + %d>)\\2 ■ |1V)I| 2 +max||n_ |^)|| 2 • ||H_ %il>)\\ 2 (142) 

m m 

< 2e (143) 

where the maximization is taken over all normalized nuclear states \<j>). We have also used Eq. (136) and the fact that 
U and II± commute. Similarly, we can obtain, 



HUfrtflMlla > v/l-46 2 . (144) 

Using Eqs. (143) and (144), we can now bound the z-polarization of the electron starting from an initial state which 
is polarized. First, we note that the operator S z can be written in terms of the projection operators 11^ and IIjj., 

S z = 5(n t -nn)® l! (145) 
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Using Eq. (144), we find that if the initial state of the system is given by ip), then the z-polarization at an arbitrary 
time t is bounded by, 



> | (1 - Se 2 ) . 



QED 



(146) 
(147) 
□ 



Proof of Theorem 5. From the statement of Theorem 5, we have defined AE max = aAE, a > 1. We have also defined 
the set of eigenstates M such that for every eigenstate k in M, we have \E k — (E)\ > AE max . Then we can write 



So, 



QED 

Proof of Theorem 6. 



(AEf) = 5> fc (25 fc -<25» s 



U\ 



keM k£M 

> J2p k (E k -(E)f 

keM 

> y^Pk (AE max f 

keM 

= (AE max f Y Pk 

keM 



... Or 



keM 



EiE k t 
Pk exp — 



> 



> 



> 



> 



EiEkt ^ iE k t 
Pk exp — ^— + 2^ Pk exp -j- 

iEut 



k^M 

EiEkt 
Pk exp — — 



keM 



E P k eX P 

keM 



h 



i{E)t ^ i{E k - (E))t 
exp — — 2^ P k ex P 1 



£ 

k<tM 



Pk exp ■ 



k$M 

iAE k t 



OL 



\ 



AE k t\ . AE k t\ 1 



> - 



> 



cos ■ 



7i a 2 



(148) 
(149) 
(150) 
(151) 

(152) 
(153) 

(154) 
□ 

(155) 
(156) 
(157) 
(158) 
(159) 

(160) 
(161) 
(162) 
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where the last inequality holds provided that, 



QED 



t < 



Hit 



2aAE 



(163) 
□ 



Proof of Theorem 7. Since } is an eigenstate of H, its evolution will correspond to multiplication by some overal 
phase exploit. Then the remainder of the proof is trivial 



U 



^li) = Ufa) 



(164) 
( 165 ) 

QED □ 

Proof of Theorem 8. Let E~ be the eigenvalue corresponding to jV^ - )- Let us consider the normalized trial state, 

\<k) = |n-,^,i) + l*.^,i> 

where 



(166) 



\h,i) 



-V 



hi ■ 



We select this trial state because it has the desirable property that 11^ |<^) 
component of \<j>) is designed to minimize the overall energy width of \(f>). We then obtain, 



(167) 



ft', V'ii ,• ) • Our selection of the spin-down 



H\<j>i) = (ffo + v)(|ft,Vfo) + _ e }_ r v 

= H^^,)-E7—l 



-E 



-E7 - Ho 
1 



-V 



ii-.i. 



V- 



1 



-Er - Ho 



-V 



ii-.i. 



-E i - H 



E i - Ho 



-V 



E- — Ho 



E, - Ho 



E- - Ho 
-V 



-V 



V- 



-E, - Ho 



-V 



-E, - Ho 



where we have used the fact that R^HoR^ = —H a . If we define the operator 



w = n t r„v 



i 



VR„ + V- 



1 



-H 



v n 



fr 



(168) 
(169) 
(170) 

(172) 
(173) 



then we can evaluate the expectation value, 

(&|#|&> = -ET + (<tH\W\4>i) 

(<t>i\H*\<fH) = (Er) 2 -2Er(4 >i \W\cl> i ) + ( ( (> i \W 2 \ ( f> i ) 

[(AH) 2 ) = (^W^^-^W^) 2 



(174) 
(175) 

(176) 

These expectation values can be bounded using the L2 norm of W (see Eq. (110) for a definition of the L2 matrix 
norm). To obtain this norm, we write W in the z basis, 



W = W dlag + W ff-diag- 

Here we define Wdiag as the diagonal portion of W in the {|z)} basis, 



z j.k 



4 42 



h 4 Al 



4 (e~ + El x + h 2 AJ2) 4 (V + Et + h 2 A k /2) 



|z> (z| 



(177) 



(178) 



6 



where j sums over nuclei for which Zj = — 1 and k sums over nuclei such that z k = +1. W //_d, as is the off-diagonal 
portion of W, 



W ff^ diag = ^2 



E 



h 4 AjA k 



h 4 A,A 



A (E7 + El z + h 2 A k /2) 4 (E7 + E# >a + h 2 A 




z') (z| 



T^AjAkiAj - A k ) 



W) <*l 



(179) 



(180) 



where z, and z' are nuclear spin configurations which differ only at nuclear positions j and k such that that Zj = 
— 1, Zfe = +1 and = +1, z' k = —1. We note that the two terms in Eq. (179) are nearly equal, leading to a significant 
cancellation in Eq. (180). This cancellation is a result of the spin echo experiment, and will be discussed in Sec. IV. 



Because h z Aj, % l A k are negligible relative to E i + Eq z , we can neglect these terms in the denominator. We also 



note that 



> 



E~ + E^ 



l|W|| 2 < 



. Then using the inequality in Eq. (112), we can bound \\W\\ 2 as 

ft 4 h 6 



~^-2,max H~ 



±2,max -rt-max 



(181) 



The second term in this equality is strictly less than the first term by a factor of e (see Eqs. (43) and (44)). So we 
can write the strict inequality, 



W\\ 2 < 



ft 4 



4 {E max + ^o.max) 



A 



2, mace 



(l + e) 



(182) 



Hence it follows from Eqs. (174) and (175) and Eq. (182) that 



H 2 \<j>}-(<f>\H\ 



< 



4 (Emax + ^o, max ) 



A 2 , m ax (1 + e) 



QED 



(183) 



□ 



Proof o/ Theorem 9. Since we know the energy width of the trial state \<fii), we can apply Theorem 6 to characterize 
its dynamics 



U\ 



= e ^- i * v ^|| 0) + A .| r -.y 



(184) 



In Eq. (184), . is an overall phase corresponding to 9 in Theorem 6. Aj is the magnitude of the residual evolution, 



and is bounded by Eq. (60). In order to characterize the evolution of 

^^l) = \<t>)-%H,i)- 

Using similar reasoning as in the proof for Theorem 7, we obtain 

U 



we then note that 



'h iPii) + A, rr.) + ( Vl - A? e^V -U)\l, 



= e 



1-A 



(185) 

(186) 
(187) 



The final element in our proof is to demonstrate the three inequalities in Eqs. (85-87). Eq. (85) follows directly from 
Theorem 6. Eq. (86) follows from the fact that we can bound the magnitude of |-f|~, 0^), 



V- 



{-Er-HoY 



< e\ 



(188) 
(1891 



Finally, Eq. (87) follows from the consideration that (from Theorem 6), (r. 



= 0. Then, 



< 
- < . 
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□ 



